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Oh! 

X ■ ABSTRACT 

M'. 

The first results on next-to-leading order QCD corrections to graviton-induced pro- 
cesses in liadron collisions in models of TeV-scale gravity are presented focusing on the 
case of dilepton pair production in pp and pp collisions. Distributions in the invariant 
mass Q, the longitudinal fraction xp, the rapidity Y and the forward-backward asym- 
metry of the lepton pair are studied. The quantitative impact of the QCD corrections 
for searches of extra dimensions at hadron colliders is investigated. It turns out that at 
the LHC (v^ = 14 TeV) the /^-factor is rather large {K = 1.6) for large invariant mass 
Q of the lepton pair, indicating the importance of accounting for these QCD correc- 
tions in the experimental search for TeV-scale gravity. At the Tevatron, the K-factor 
does not substantially deviate from the Standard Model value. However, its inclusion 
is necessitated to make the cross-section stable with respect to scale variations. 

PACS numbers: ll.lO.Kk, 12.10.-g, 12.38.-t, 13.60.Hb. 
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1 Introduction 



In the recent past, models of extra dimensions have been studied as serious contenders 
for physics beyond the Standard ModeL These extra-dimensional models are attractive 
because they help address the hierarchy between the electroweak and the Planck scales. 
In the model of Arkani-Hamcd, Dimopoulos and Dvali (the ADD model) [1] the Planck 
scale is reduced to a TeV and the hierarchy problem thereby avoided by invoking a 
large magnitude of d extra dimensions. 

In the ADD model, the Standard Model (SM) fields are constrained to a 3-brane, 
while gravitons propagate in the A+d dimensions. Then the size of the extra dimensions 
is only constrained by the length scales to which the gravitational inverse square law 
has been experimentally tested, which are currently probing the sub-millimetre range. 
For d between 2 and 6, the size of the extra dimensions varies from a millimetre to a 
Fermi. The relation between the 4-dimensional Planck scale Mp and the scale Ms in 
(4 -I- d)-dimension is 

Mp ^ M^^^^^ R'^ , (1.1) 

where it! is the compactification radius. In the ADD model, because it! is large it is 
possible to lower Ms down to a TeV and avoid the hierarchy problem. An important 
consequence of the lowering of the Planck scale is that quantum gravity effects could 
be tested at energies of 0(TeV). So it opens up a new area of studies on the effects of 
gravitons at present and future colliders. 

The A + d dimensional graviton manifests in 4-dimensions as a tower of massive 
Kaluza-Klein (KK) modes. These KK modes interact with the SM particles confined 
to the 3-brane via the energy-momentum tensor. Each KK mode couples to the SM 
particles with a coupling of the order of 1/Mp, (we will use k, in the rest of the paper). 
As can be seen the coupling of each KK mode to the SM particles are highly suppressed. 
But the effective coupling is obtained after summing over all the KK modes and due to 
the high multiplicity of the KK modes the effective interaction has a strength of I/M5 
[2, 3]. This enhanced coupling provides viable signatures of the graviton KK modes at 
colliders. Both real graviton production and the effects of virtual gravitons in various 
processes have been studied in the literature [4] and have yielded bounds on Ms which 
are in the ball-park of a TeV. 

These collider processes involving gravitons have, however, been studied at the 
leading order (LO) in QCD ^. At hadron colliders such as present Tevatron and up- 
coming Large Hadron Collider (LHC), the theoretical uncertainties coming from the 
QCD effects due to initial state partons are quite sizable. The sources of these uncer- 
tainties are two fold: renormalisation and factorisation scales and the parton distribu- 
tion functions. The scale uncertainties can come from the strong coupling constant as 
well as parton distribution functions. Experience with next to leading order (NLO) 

^Recently NLO-QCD corrections to e+e~ — > hadrons have been studied in [5] 
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contributions to SM processes strongly suggests that the LO corrections at hadronic 
colliders are quiet unreliable [6, 7]. The standard Drell-Yan (DY) processes and Higgs 
production at hadron colliders, for example, not only get large corrections but the 
theoretical uncertainties get reduced significantly from NLO corrections. 

It is with this motivation that the present paper presents the results of the com- 
putation of NLO-QCD corrections to the DY process, Pi + P2 — /^"""A*" where Pi, 
P2 are initial hadrons and X is an arbitrary hadronic final state. This process takes 
place via the exchange of 7, Z and graviton. These NLO results are used to study 
the impact of the QCD correction for this process at the present Tevatron collider and 
the future LHC which is planned to operate at centre-of-mass energies around 14 TeV. 
This process had been considered earher to LO [8, 9] in the ADD model. 

The results we present here are for the ADD model, but since the QCD corrections, 
presented here, are model-independent they may equally be used for studying the 
Randall-Sundrum model of warped compactification [10]. 

This paper is organised as follows. In Sec. 2 we describe additional graviton vertices 
needed to perform the NLO QCD corrections. We study the distributions in the 
invariant mass Q (Sec. 3), the longitudinal fraction xp (Sec. 4), the rapidity Y (Sec. 5) 
and the forward-backward asymmetry (Sec. 6) of Icpton pair. Finally in Sec. 7 we 
present the discussion and summarise the results. In Appendix A and B we present the 
detailed expressions needed for xp and Y distributions in Sec. 4 and Sec. 5 respectively. 
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2 The Model 



We work with the following action 



(2.1) 



where Ssm is the Standard Model action, h^" is the graviton field and k is the strength 
of the interaction. The energy momentum tensor in QCD is given by 

+^{Atd,{d-A-^) + A^d^id'^A-^)) + - igT'^AD^p 



(2.2) 



In the above equation, ^ is the gauge fixing parameter. We work in the Feynman gauge 
in which the gauge parameter ^ = 1. We have displayed explicitly the ghost terms with 
the ghost fields u;°'{x) since they contribute to our one-loop virtual corrections to the 
process under study. The presence of the ghost fields introduces two new vertices (see 
Fig. 1). 

1) graviton- ghost-ghost vertex 



G 




C, P:P3 



a, Pi 



Figure 1: Ghost vertices 



(2.3) 
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2) graviton- ghost-ghost- gluon vertex 




(2.4) 



where 



(2.5) 



The indices refer to the graviton. The momenta pi,p2 and the colour indices a, h 
correspond to the incoming and outgoing ghosts respectively. Finally c, p, indicate 
the colour index, Lorentz index and momentum of the gluon. The other Feynman rules 
are given in Ref. [3]. ^ 



^The only exceptions are the Feynman rules for the fermion-anti-fermion-gauge boson-graviton 
vertex and the three gauge boson-graviton vertex which differ from those of Ref. [3] by an overall 
sign. 
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3 Invariant lepton pair mass distribution da jdQ^ 

We start by considering Pi, P2 scattering to leptonic final states, say /i"*", [x~ 

PliPl) + P2iP2) ^ l^^ik) + l^'ik) + XiPx) , (3.1) 

where pi,p2 are the momenta of incoming hadrons Pi and P2 respectively and p~,fi~^ 
are the outgoing leptons which have the momenta /i,/2- The final inclusive hadronic 
state is denoted by X and carries the momentum Px- In the QCD improved parton 
model, the hadronic cross section can be expressed in terms of partonic cross sections 
convoluted with appropriate parton distribution functions as follows 

2^^^(T,g') = E J dx,j dx2 f:\x,) f^\x2) 

V ab=q,q,g ^ ^ 

pi J^ab 

X dz 2s [z, 6{t - ZX1X2) . (3.2) 

The scaling variables are defined by ki — Xipi, k2 — X2P2 where ki, k2 are the momenta 
of incoming partons. 

(P1+P2)' = S, {k, + k2)^ = s, {h + hf = q.q = Q\ 

T = — , , T^XiX2Z. (3.3) 

o s 

m 

The partonic cross section for the process a{ki) + 6(/c2) j{~<l) + '^Xi{~Pi) is given 

i 

by 

da"^ 



In the above equation, the sum over Lorentz indices between matrix element squared 
and the propagators is implicit through a symbol "dot product" . The m + 1 body 
phase space is defined as 

m 

x(27r)"5(")(A;i + ^2 + g + EP^)' (3-5) 
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where n is the space-time dimension. The propagators are 



Pg(?) = V{Q')B,,^,{q) = B^,^,{q)PG{Q') , (3.7) 



where 



2 

Th -L 

= -9,^u + ^ . (3.8) 
The summation of the virtual KK modes in the time-hke propagators [3] leads to 

where the integral I is regulated by an ultraviolet cutoff, presumably of the order of 
Ms [2, 3], which sets the limit on the applicability of the effective theory. For the DY 
case under consideration this consistency would imply Q < Ms- Further relating the 
gravitational coupling, the volume of extra dimension and the cutoff scale [3] 

K^R'^ = 87r{47r)'^/^T{d/2)Ms^'^+^^ , (3.10) 

we express the function T>{Q'^) as 

The summation over the non-resonant KK modes yields 

d/2-l . . 

H^) = - E - 2 - 1) ' d = even, (3.12) 

k=l 



m = - 'T S^-^'-' + \ log (^) . = odd . (3.13) 

K=l 
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where uj — Mg/Q- The leptonic tensor involves the computation of square of the 
matrix element for the process ^/Z/G ^ /i'^ + /i" and the 2-body phase space integral 

C"'-^^^-{q) = / n ((C|^27r5+(/,2)) (27r)"5(")(g -h- k)\M^^'-^^'^- f . 

(3.14) 

The leptonic part is easy to compute. It is equal to 

£.r->M+M-(^) = ^^,{q)L,AQ') , jf = {^^,ZZ,^Z}, 

CGG^,^,-^q) = B,.,M)Lg{Q^). (3.15) 

with 

where a is the fine structure constant 

Cw — cos 9w , Su, = sin , 

^ = ^ra^-4Qa, ^f = 4^a^ (3.17) 

and Qa is the electric charge of quarks and leptons. Hence, we obtain 

^-^^Tn^^^'^') = ^ ^ P;(g^)i^^(g^)L,,,(g2)H^^!^^(r,g2). (3.18) 

The hadronic structure function is equal to 

W,5^^(r,g^) = j dx, j dx, f:^{x,) f^^{x,) 

X Jdz 5{r - ZX1X2) J dPS^+i \M^^^"'\^ Tjy{q) , (3.19) 

with 

Tjj' = Vn'^iQ), jf ^ll,lZ,ZZ , 

Tgg = B^,p4q). (3.20) 
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Figure 2: Born contributions: 

The tensors fJnuiq) and Bf^ypfj{q) are defined in Eq. (3.8). To compute the distri- 
bution of the di-lcpton pairs, the matrix element squared \M°''^^^^ p Tjji{q) has to be 
substituted in Eq.(3.19) provided the integrations over dPSm+i and dz are performed 
in a suitable frame. We define the bare partonic coefficient function A^j, {z, Q"^) as 

A^iV, Q') = j dPSrr.^, \M'^^-^^'\' Tjjiq) , (3.21) 

where 

g2^2 -^^ V— ; 

There are two classes of processes that contribute to the partonic cross section. The 
first one has only a virtual photon or a Z boson whereas the second one only contains 
a graviton in the intermediate state. Interestingly, the interference between first and 
second class diagrams identically vanish when the phase space integration is performed. 
In the case of the photon exchange process, it is easy to understand if we realise that 
there is no third rank tensor, say S^'^f which can be constructed out of gi^^ and Qp 
(q.q 7^ 0) (which are the only tensor and vector available in our problem) that satisfies 
qp,S^'^^ = (as it should be for the theory where we have gravity coupled to a conserved 
energy momentum tensor). We found a similar argument for the Levi-Civita tensor 
which shows up in the case of the electro- weak vertices. Therefore also here there are 
no interference terms. We could also verify this by explicit computation. At LO the 
following process contribute when a photon or Z-boson appears in the intermediate 
state 

q + q^ 1*/Z* . (3.23) 

In the case the Graviton appears in the intermediate state, we have two processes at 
LO (see Fig. 2) 

q + q^G*, g + g^G*. (3.24) 



jf = ZZ, -fZ, 
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In NLO for the photon or Z-boson we have the following 2 — > 2 reactions [6, 7, 11] 

q + q ^ /Z* + q, g + g ^ 7*/^* + one loop, 

q + g^^*/Z* + q, q + g^Y/Z* + q, (3.25) 
and for the graviton we have (see Figs. 3-9) 

q + q G + q , q + q ^ G + one loop , 



q + g^G + q, q + g ^ G + q, 

g + g G + g , g + g ^ G + one loop . 



(3.26) 



The cross sections beyond leading order involve the computation of one loop virtual 
gluon corrections and real gluon bremsstrahlung contributions to leading order pro- 
cesses. We also include processes with a gluon in the initial state. Since we are deahng 
with energy momentum tensor coupled to gravity expressed in terms of renormalised 
fields and masses, there is no overall ultraviolet renormalisation required. In other 
words, the operator renormalisation constant for the energy momentum operator is 
identical to unity to all orders in perturbation theory. But we encounter soft and 
coUinear divergences in our computation. We have used n-dimensional regularisation 
to regulate both these divergences. To that order, we have defined n — 4: + e where 
n is the space-time dimension. With this procedure all divergences appear as 1/e" 
where a = 1,2. The soft divergences coming from virtual gluons and bremsstrahlung 
contributions cancel exactly according to the Bloch-Nordsieck theorem. The remaining 
coUinear divergences are removed by mass factorisation which in our paper is performed 
in MS scheme. The Drell-Yan coefficient function after mass factorisation, denoted 
by Kbi^^ <5^ /^f)' is computed by 

Ai,{z, Q\ l/e) = V^{z, fil, 1/e) r,,{z, fi^, 1/e) ® Ai,{z, Q\ ix\) , (3.27) 

where A^^(z, Q^^l/e) is the bare partonic coefficient function before mass factorisation 
is carried out. Further we have dropped the double index jj' because of the vanishing 





G 



Figure 3: Virtual corrections 
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Figure 4: Virtual corrections, gluon loops. 

9 

G 




9 



G 




G 



Figure 5: Virtual corrections, quark loops. 



interference terms and replace it by the single index i instead. The factorisation scale 
is given hy and ® is the convolution symbol defined as 



and the kernel Tcd{z) in the MS scheme is given by 

e 



(3.28) 



(3.29) 



9 



G 



G 
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Figure 6: Virtual corrections, ghost loops. 
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Figure 7: Real emission, q q ^ g G. 



where Pcd\^) the leading order Altarelli-Parisi sphtting functions [12], 



A » _ A (0),i , A (l),i 



For convenience we define as 

Expanding Eq. (3.27) up to order a^, we find 



Ajf = a4r£)»A(r^ + a.A(W-, 

A« = Ajf + a,HrW»ASf + a.Ar, 

A° = a,i(C»Ar + r«»A<r)+a,A<i)«. 



(3.30) 



(3.31) 



G 



(3.32) 
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G 




G 




G 




Figure 8: Real emission, q g ^ q G. 

Prom the above expressions one can compute the coefficient function lS,\^{z,Q'^ , Hf) 

from the bare ~S'g^{z,Q'^ , jip, l/^) and the known AltareUi-Parisi kernels . Finally 
wc have to fold these finite /S.\{z,Q'^ , jip) with the appropriate partonic distribution 
functions to arrive at the distribution for the DY pair. For completeness we present 
the results below 



da^^^'^ /•! rl rl 

2S ——^{t,Q'^) = ^^SM,qJ^ dxi dx2 dz S{t - ZX1X2) 



dQ 



X 



H,,{x^, X2, 111) (a(J'^/^(^, Q\ 111) + a,A«'^/^(^, Q\ iil)) 



+ [Hqg(xi,X2, Hl) + Hgq(xi, X2, //|.)) A ^'^/^ 



+ / dxi / dx2 / dz 6{ 

„ Jo Jo Jo 



T - ZX1X2 



X 



Hgg{xi, X2, 111) (aJJ'^(z, Q\ 111) + a,Ag'^(z, Q\ iil)) 
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9 




mm 




Figure 9: Real emission, g g ^ g G. 



-rEgg{x^, X,, 111) (a£)'«(;., g^ + a,AW'«(^, g^ /x^) 



(3.33) 



The constants TsM,q-i are given by 
4a2 



SM,q 



3Q2 



G 



2m2 



3207r2 ' 

The renormahsed partonic distributions are 

Hqq{xi,X2, 1^%) = fq'{Xl, n]r) fp{x2, 1^1) + fg'{Xl, 1^%) fq^{x2, fll) , 



(3.34) 

(3.35) 
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Hgg{Xl,X2,l4p) = Hgg{X2,Xl,IIp) , 

Hgg{Xl,X2,^ll) = fg^{Xl,^l) fg^{X2,^il) . 

The LO coefficient functions are 



'■qq 



2n 

TV 



S{l-z), 



(3.36) 



,(0),G 
^qq 



TT 



a(o),g 

99 



and the NLO results are given by 



TT 



2{m - 1) 



,{1)7/2 
^99 



(_4 + 2C(2))5(l-.)+2j^ln(|) 



(3.37) 



+ 



'l-z 



^q{q)9 



^qq 



2{l - 2z + 2z')\yi 



Hlz 



l + 6z- 7z' 



+4 



/ ln(l -z) 



-2—^ H 

1 — z 3z 3 



a(1)G 



2(-7+g + 2. + 2.-)lnf^^V-)" 
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12 , 

+9 + lSz- 

z 



99 



^)c.[(-H5!+8C(2)y(l-.)+8^1n(|) 



2(iV2 



[-2 + - + Z- z^) 
z 



X In 



+ 



TT 



ln(2) „ 22 „ 22^2 



^2(Ar2 - 1) 

For SU (N) the colour factors in the above equations are 



(3.38) 



Cp 



2N 



Ca^N, 



(3.39) 



and n/ is the number of flavours. The "plus" functions appearing in the above results 
are the distributions which satisfy the following equation 



'\n\l - z) 
1-z 



'ln*(l-z)^ 



1-z 



£dzU{z)g{z) = £dzf{z){g{z)-g{l)), 
where g{z) is any well behaved function in the region < ^ < 1. 



(3.40) 
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4 Differential cross section with respect to xf 

In this section we compute the xp-distribution of the di-muon pair up to NLO in QCD 
with the 7, Z and the graviton in the intermediate state. The NLO corrections to 
processes containing an intermediate photon were already calculated in [6], [13]-[16]. 
The variable xp is defined as 

xf= ^ , (4.1) 

where Pi,P2 are the momenta of incoming hadrons and q is the sum of final state muon 
momenta, q — li + l2- In the CM frame of the hadrons (CMH) 

XF - , (4.2) 

where qcuH is the third component of q in the CM frame of Hadrons. The hadronic 
xp distribution can be obtained by introducing the identity 

dx,6 (x, - ^^i^^Htl^ = 1 , (4.3) 
in Eq.(3.2) and bringing the measure dxp to the left. Hence we find 

^ jf=j,Z,G ^ 

(4.4) 

The hadronic structure function is defined as 
dW^'f' 



{t,xf,Q') = E /^^i/ dx2 f^Kxi) f^\x2) 

ab,jj' 

X J dz J dPSm+i |M"^^^'^'|2 Tjfiq) 6{t-zxiX2) 

^S^,,-^JPlzMliy (4.5) 

We start by computing the Born contribution. This involves the computation of the 
matrix element squared for process a(/ci) + b{k2) — > jj'{—q) and the (0 + l)-body 
phase space integration as well as the dz integration. The (0 + l)-body phase space is 
computed using Eq. (3.5) and one finds 

J dPSo+i = ^5{z - 1) . (4.6) 
17 



To compute the Born process a{ki) + b{k2) — > jj'{—q), we choose the following param- 
eterisation for the momenta 

A;i = ^(l,0,---,0,1), A;2 = ^(1,0, •••,0,-1), -g = VI(1, 0, • • •, 0, 0) . (4.7) 

Subsequently we express the partonic momenta in terms of the hadronic momenta 
using the scaling variables xi,X2 as ki — xipi, k2 — X2P2 

J dPSo+i J dz5 (^F - ^i^L_^H]_£j _ ZX1X2) = ^ J dz5{xF - xi+ X2) 

x5{z — 1)S{t — ZX1X2) ■ 

(4.8) 

We choose to work with the variables x?, X2 which are defined through 

XF=x'i-x'^2, T = x[xl. (4.9) 

Solving the above equations, we find 

4 = li^F + \lxl + At) , xl = ]^{-XF + \lxl + At) . (4.10) 

After rewriting the arguments of the delta functions in terms of the variables x\,X2-, 
and performing the z integration using 5{t — ZX1X2), we arrive at a simple looking 
expression for Eq. (4.8) as 

JdzSiz-mxF-x. + X2)SiT-zx.X2)\M^^-^^r-T,^- 



Q 



27^ 6{x,-x^,)6{x2-xl) ^^j^, 2 , . . 



Finally, the Born matrix element squared 

|M"^-"^'^"|2 •% , (4.12) 

is computed using the parameterisations given in Eq. (4.7). It is then substituted in 
Eq. (4.5) to arrive at the leading order xf distribution. Our next task is to compute 
the next to leading order contributions to the Born xf distribution. This involves 
the computation of the matrix element squared for the processes a{ki) + b{k2) — > 
jj'{—q) + c{—k) and 1 + 1 phase space integration. Since the integral on the right hand 
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side of Eq. (4.4) is Lorentz invariant, we can work in the CM frame of the incoming 
partons. In this frame the momenta of the particles are parametrised as 



h = ^(1,0, •••,0,1), A;2 = ^(1,0, •••,0,-1), 



V S 

—q — ^{1 + z,0,- ■ ■,—{1 — z)sin9,—{l — z)cos9) , 



-k = ^{1- z,0,- ■ -,{1- z)sme,{l- z)cos9) . (4.13) 
In this frame, the 1 + 1 phase space becomes 

/ "^^^ = h Tj^r''^' - * - f • ("^) 

where y is related to cos^ by cos^ — 2y — l. The xp delta function becomes 

5 (^^ _ 2(pi -V2) ■ g j ^ ^ _ ^ ^ _ ^ _ (4_;^5) 

Wc can perform y as well as z integrations by writing the delta functions in terms of 
the variables 

5{x,-x, + zx, + y{l-z){x^ + x,)) = + "^'^ _ ,o,0) ^ " ^1 ^ 

5(T-^a;ia;2) = ^— 5(z-z*), (4.16) 
X1X2 



where 

^* ^ X2(xi — xf}{xi + x") ^ ^l3^2 
{xiX2 - Xlx^){xi + X2) ' X1X2 ' 

In other words 



I dPS,+, I dzS{r - zx,X2)5 ^XF - /'^"^ ) \M''''^"'\^ ■ Tjj, = 

l_ (Q^Y^' 1 (a;;x°)-i(xi + X2)-^ 
87rl^47r; r{l + e/2) xiX2ixi + X2) 

X ((xi - x',){x^ + a:°)(a:2 - 4)(^2 + x[)f \M'^'-^^'\' ■ 



y=y*,z=z* ■ 

(4.18) 
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Notice that the matrix elements \M"'^^^^'\'^ • Tjjt are evaluated for y and z given by 
the delta function constraints given in Eq. (4.16) and Eq. (4.17). To this order, one 
has to include the virtual corrections to the Born processes as well. We use the same 
formula given in Eq. (4.6) with the one loop corrected matrix elements. The sum of 
real emission contributions and the one loop corrections to the Born processes is free 
of the soft singularities as expected. Therefore the result contains only the coUinear 
divergences which can be removed by the standard mass factorisation. This is done by 
redefining the parton densities using the Altarelli-Parisi kernels as follows 

/<r(^)=Er„l.'®/r(^>/^F), (4.19) 

6 

which implies 

/;(^) = /;(^, 1^) - <^s-nf (r« ® f^{z, 111) + ® f!{z, 111 

-<^s-r(^J®f^i^,^l)- (4-20) 



Finally we arrive at 



^ SM,q ( Cg-^(Xi, X2, III) + Cq^(Xi, X2, nl) + Cg^{x\, IlI) 
i=q \ J 

+ E -^G (c%{x\, x% III) + Cg{xl xl 111) + C^^ixl xl 14)] 



i=q 



+J^GC%{x'„xliil), (4.21) 



where 



^abi^l^ ^2) A*f) — ^ab^ ^2) A*f) + C'afe'^ ^2) A*f) ■ (4.22) 

We have presented C^^'^'^\ C^^'^^\ C^'*-""* and C^'*-"^'' in appendix A. 



20 



5 Differential cross section with respect to Y 

In this section we compute the rapidity distribution of di-lepton pair up to NLO in QCD 
with 7, Z and graviton in the intermediate state. The NLO corrections to processes 
containing an intermediate photon were already calculated in [6], [13]-[16]. We define 
Y as 

2 \Pi-qJ 

In the CM frame of the hadron (CMH), this is nothing but 

Y ^llog (^^MKZ^^Y (5.2) 

2 V QCMH + QCMH J 

With this definition the distribution in Y can be computed by introducing the identity 

in Eq. (3.2) and bringing the measure dY to the left. Hence we find 

25 1^(^^Y,Q^) = 2^ ^ ^^-(^'^ ^^-^'^^'^ ^l^^^'^'^'^ -(^-^^ 

^ n'=i,z,G 

The hadronic structure function equals 

dwr^''' 



X J dz J dPSm+1 \M''^^"'\^ Tjj,{q) 6{t-zxiX2) 

We first compute the matrix element squared for a + 6 — > jj' and then substitute this 
in Eq. (5.4). Then we perform the (0 + l)-body phase space integration as well as z 
integration. We follow the same steps as for the xp distribution. Using Eq. (4.6) for 
the (0 + l)-body phase space and performing z integration, we obtain 

/dPSo. ji^s [y - ilog 5(r - ^ (y - llog (I)) 

x5{z - l)5(r - ZX1X2) . (5.6) 
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Instead of working with the variables Y,t, we choose to work with x1,X2 which are 
defined as 

^ = ^log(|), r^xlxl (5.7) 

Solving the above equations, we obtain 

xl = V^e^ , 4 = v^e-^ . (5.8) 

We can perform the z integration in Eq. (5.6) using S{t — ZX1X2) after rewriting the 
remaining delta functions in terms of the variables x?, x^. We then get 

^ / dz6{z - 1)S (y - 1 log (^)) 6{r - zx,x,)\M-'-"' \' = 

^5(a;i - x\)6{x2 - xl)\M'^'-^^' \' T,y|,=i . (5.9) 

Next we evaluate the NLO correction to the Born processes. This involves the compu- 
tation of 1 + 1 phase space for the process a{ki) + h{k2) + c{—k). Since the 
integral on the right hand side of Eq. (5.5) is Lorentz invariant, we can work in the 
CM frame of the incoming partons. We follow again the derivation of xp distribution. 
The Y delta function in this frame becomes 

The two delta functions reduce to 



1 f xi{l — y{l — z))\\ 2x1X2x5x2(^1X2 + x^xl 



0^0-) 



2 \X2{z + y{l - z)) J J (xiX2 - X?X^)(xiX^ + X?X2) 

5(t — ZX1X2) — S (z — z*) , (5-11) 

a:;ia;2 



where 

(X1X2 + X2Xi)(xiX2 - X1X2) ' X1X2 ' 

Using Eqs. (5.11) and (5.12) we arrive at 



J dPS,+, J dz 5{t - zxyX2) ^{^ - \ log (^)) l^"'"''T T^j' = 

1 ( ^ 2x^X2 (XiX2 + x'j'x2) , _|_ \-e 

8^\I^) r (1 + e/2) xiX2(xixO + X2X?)2 ^^'^^ + ^1X2) 

((xi - x?)(xi + x2)(x2 - x°)(x2 + x°))'^' IM'^'-")' T^Ay=r,z=z* (5.13) 



22 



We substitute Eq. (5.13) in Eq. (5.5) to obtain the real emission contributions to the 

Y distribution. The virtual corrections to the Born processes can be obtained using 
Eq. (4.6) with the one loop corrected matrix elements. After adding the real emission 
contributions and virtual corrections to the Born processes, all soft singularities can- 
cel. The remaining coUinear divergences are removed by mass factorisation using Eq. 
(4.20). Finally we arrive at the finite results for the Y distribution 

+TgD%{xIxIiiI), (5.14) 

where 

-^a6^(^l' ^2) A*f) — ^ah ^ (^1) ^2) A*f) + \x\, X^, 1^^) , 

Dg{xlxlf,l) = D^f\xlx',,pil) + a,D^^^'\xlxlpil). (5.15) 
We have presented D^^'^^\ D^^'^'^'^^ -0^/°'* and -D^j,'*^^^ in appendix B. 
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6 Forward Backward asymmetry Aps 

In this section, we compute the forward backward asymmetry of the di-lepton pair. 
The NLO correction to ^4^^ in the SM were calculated in [19]. The forward backward 
asymmetry is computed using the distribution define by 

(-■«^) - (i -/_ j 1^ (^'^ ■ 

where cos^^* is the angle between the final state lepton (say li) momentum and the 
initial state hadron (say pi) momentum in the CM frame of the lepton pair. In general 
it can be written as 

Pi ■ {h - h) ,„ „^ 

cosfc* = , (6.2) 

The computation of the above distribution is quite similar to that of da/dQ^. In the 
QCD improved parton model, we find 

^ ai>=q,q,g 

^-^28-— -(z,Q^)5{t-zx^X2). (6.3) 

Jo dQ^dcosO* ^ ' 



We compute the partonic cross section using 
dn°^ 1 /• 

im^' - Y. .^.J ■ ^^i'* ■ ""M 



jr,jj'^^^i-{q^ cos 9*), (6.4) 
where O^'^^^^ i^q^ cos 6**) can be computed using 



f=iV(27r)" 

xs (cos 0* - \M^''^^^^~ r ■ 

V Pl-(/l+/2)/' 



(6.5) 



The computation of the leptonic part using the above formula is straight forward. The 
hadronic part again involves the computation of various processes that contribute to 
Q distribution. Since we are looking at the angular distributions which are "odd" in 
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cos^*, the contributions come mainly from interferences. The non- vanishing contri- 
bution in the standard model sector comes from the interference of photon mediated 
processes with Z-boson mediated processes. We also find that non-vanishing con- 
tributions come from the interference of standard model diagrams with the graviton 
exchange diagrams. These inference diagrams are absent in the computation of Q, Xp 
and rapidity distributions because they are odd in cos^*. We have regularised all the 
divergences using dimensional regularisation. The remaining mass singularities are 
removed by the mass factorisation. To the end, we find 



y^^SJ^SM,q f dXi [ dX2 [ dz 6(t - ZX1X2) 

n -fO Jo Jo 



X 



+8Hgq{xi,X2, jjil) ias6^^g^'^^{z, nj. 



+ 



SJ-'a / dxi / dx2 / dz d{T — ZX1X2) 
„ Jo Jo Jo 



X 



SHgg(xi,X2, i^f) [6A'g^'^{z, /x^,) + Q ^5 A^^' '^^ {z , Q'\ jip 



(1),G/ ^2 ,,2 



+5Hgg{xi, X2, i^f) (asSA''gJ''^{z, 



+SHgg{xi, X2, nD) (a,5AW'^(2;, 



The constants SJ-'sM,q, are given by 

(Q2 _ Ml) 



5Tt 



SM,q 



QqQegta^ 



A A 
q He 



(6.6) 



■((g2-M|)2 + M|r|)c44^^^^^«^^ 



47r 



5Ta = 



QqQe ~l~ 



Q\Q'-Mi) 



((g2_M|)2 + M|r|)4.2 



V V 

9q 9e 



(6.7) 
(6.8) 
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The renormahsed incoming partonic fluxes are defined by 



5Hgq(xi, X2, = fg'(xi, ^p) (/f'(x2, - fp{x2, nl)) , 

SHqg{xi,X2,f4') = SHgq{x2,Xi,II%) . (6.9) 

The LO coefficient functions corresponding to Eq. (6.6) are 



SA 



2n 
TV 



TT 



The NLO contributions are given by 



5A 



(1)7Z 
99 



^99 + TV 



4(1 + ^) ln(^)+4(l-z) 



(6.10) 



99 



99 



AT- 



4(1 - z) \ii{z) 



m 



8^^ In(^) + 2(1 + 2z - 3^2) 



5A, 



;i)G 

99 



1 a{1)G , 



4(1 + ;2)ln(;2) -25(1- ;2) 



"^99 16^99 ^ f 



4(l-2z)ln(^)+4(l-z) 



5A(i)G 

99 



.1a(^)^ + — 
16 8Ar 



8(2 - ;2 + ^2) in(^) + 2(1 + 2z- ?,z^) 



. (6.ii; 



The coefficient functions A^^^'* with i — ^,Z,G correspond to the ones found in the 
cross section for the invariant lepton pair distribution Eq. (3.33). They can be found 
in Eq. (3.38). The forward backward asymmetry is defined by 



2 9 (r 



(6.12) 



The impact of the NLO correction is discussed in the next section. 
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7 Discussions 



In this section, the effect of the NLO QCD corrections on various distributions such 
as the invariant mass Q, the longitudinal moment fraction xp and the rapidity Y of 
the lepton pair are presented. We have chosen to present these distributions for the 
LHC (v^ = 14 TeV) and for the Run II of Tevatron (v^ = 1.96 TcV). Wc want to 
emphasise that we are not analysing the existing Tevatron data to extract bounds on 
Ms, which will require doing a full hadron-level simulation, but are only interested in 
gauging the impact of the QCD corrections on the bound. 

The standard model parameters which enter our analysis are a — 1/137.03604, 
Mz = 91.1876 GeV, Fz = 2.4952 GeV and sin^ Ow = 0.227. For the ADD model 
parameters, we choose Ms = 1-5(2) TeV and d = 2(3) for the Tevatron (LHC) re- 
spectively. For the parton density sets we adopt in leading order MRST 2001 LO 
(A = 0.1670 GeV) and in next-to-leading order the MRST 2001 NLO (A = 0.2390 
GeV) [17]. For LHC and the Tevatron we choose the range 150 GeV < Q < 1100 GeV. 
Further wc have for the LHC the ranges —0.7 < xp < 0.7 and —2.7 < Y < 2.7 at 
Q = 700 GeV. A similar analysis is performed for the Tevatron. Here wc choose for 
the xp and Y distributions the regions —0.7 < xp < 0.7 and —0.9 < Y < 0.9 at 
Q = 400 GeV respectively. The renormalisation scale and factorisation scale np 
are taken to be equal. Finally the factorisation scale /ip is chosen to he /ip — Q unless 
mentioned otherwise. 



Let us first study the invariant lepton pair mass distributions 

da\Q) 



(7.1) 



dQ 

ioT I — SM (standard model) and / = GR (gravity) and the sum of both mechanisms 
SM + GR. 

Starting our discussion with the LHC we have plotted the distribution in Eq. (7.1) 
for 150 GeV < Q < 1100 GeV in Fig. 10a. The standard model dominates the cross 
section for Q < 600 GeV but for Q > 600 GeV the graviton mediated process takes over 
and overwhelms the standard model result. The dominance of the graviton becomes 
very clear in the longitudinal momentum fraction cross section Fig. 10b 

da\Q, Xp 



dQ dxp 

and the rapidity distribution in Fig. 10c. 

da\Q,Y) 



dQ dY 



(7.2) 

Q=700 GeV 



(7.3) 

Q=700 GeV 



In a large invariant mass bin like Q = 700 GeV, the effect of the graviton is clearly 
discernible in both the xp and Y distributions especially for the central xp and Y 
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regions. At large absolute values for both variables the cross section for the graviton 
mediated processes is again smaller than the standard model result. 

In Fig. 11a we have studied the K- factor for the invariant lepton pair mass distri- 
bution defined by 



\dalo{Q)] 


-1 


d(7NLo{Q) 


dQ 




dQ 



where NLO stands for the next-to-leading order corrected cross section. A similar 
definitions exist for the longitudinal momentum fraction in Fig. lib 



daioiQ.xp) 



n -1 r 



dQ dxp 
and the rapidity in Fig. 11c 



{Q,Xf) 



dQ dxi 



(7.5) 



Q=700 GeV 



dQ dY 



n -1 r 



d a^^^o{Q,Y) 
dQ dY 



(7.6) 



Q=700 GeV 



Such a definition is useful for Q, xp and Y distributions because there is no interference 
between the standard model and the graviton mediated processes. We find that is 
moderate for all values of Q in Fig. 11a. On the other hand K*^^ is much larger than the 
standard model value for K^^ . In particular for Q > 700 GeV, K'^^ overwhelms the 
standard model result K^'^ completely so that the whole X'S'a^+G-R (.g^^ ]-,g attributed 
to K*^^. This is because the graviton mediated processes show up already at the 
Born level (see Eq. (3.35)) and at the LHC the gluon flux is quite large which lead 
in NLO to big effects mainly due to small x terms in the coefficient function A^^^*^ 
(Eq. (3.38)). These effects are also visible in Fig. lib for the xp and in Fig. 11c for 
the y -distributions. If the choice of parameters Ms — 2 TeV and d = 3 is reasonable 
then for Q > 700 GeV the effect of the graviton will be observed in view of the cross 
section which amounts to about 10^^ pico barn. 

In Fig. 12a, Fig. 12b, and Fig. 12c, we have plotted the scale variations of the 
various distributions for both LO and NLO cross sections. For this we define for 
the invariant lepton mass distribution in Fig. 12a as 



Rlo — 



daioiO^ji^^io) 
dQ 



1 -1 



dQ 



(7.7) 



Q=700 GeV 



r{ 



NLO 



dQ 



d(^NLo{Q^^) 

dQ 



(7.8) 



Q=700 GeV 
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where fio is a fixed scale which is chosen to he /jLo — Q. For the longitudinal fraction 
of the lepton pair we have in Fig. 12b 



Ri 



LO 



daio{Q,xp,fi = /ip) 
dQ dxp 



daio{Q,XF,fi) 
dQ dxp 



Q=700 GeV,a;F=0 



(7.9) 



^NLO 



d^NLoiQ-n-- fl = /io) 



dQ dxi 



d(JNL(){Q-''t-J^) 

dQ dxp 



(7.10) 



Q=700 GeV,XF=(i 

Finally we have plotted the scale variation for the rapidity in Fig. 12c 



dcrio{Q,Y,ii^ IIq) 
dQdY 



d(Tlo{Q,l^) 
dQdY 



Q=700 GeV,Y=0 



R 



NLO 



dQdY 



n -1 r 



d(TNLgiQ^ 
dQdY 



(7.11) 



(7.12) 



Q=700 GeV,y=0 



In Fig. 12a we observe that for the invariant lepton mass distribution the scale variation 
is appreciably reduced for < fiQ while going from LO to NLO. For /x > /xq the same 
is happening but here the effect is not very big. However the inclusion of the standard 
model into gravity makes the situation better in particular for fi < iiq. This holds 
for LO as well for NLO. Apparently the gluon initiated processes in the case of the 
gravity show a larger scale variation than occurs in the case of the standard model. 
If we look at Eq. (7.9,7.10) plotted in Fig. 12b the features are the same. There is 
an improvement in scale variation while going form LO to NLO. The same holds for 
Eq. (7.11, 7.12) plotted in Fig. 12c. 

Now we turn our attention to the Tevatron. In Fig. 13a we show the result for 
the invariant mass distribution of the lepton pair. Here the graviton mediated process 
becomes larger than the standard model result at Q > 700 GcV. In the xp-distributions 
(Fig. 13b) and the F-distributions (Fig. 13c) the gravity effects are small. The picture 
for the i^T-factor differs from the one observed at the LHC In Fig. 14a the i^T-factor 
for the graviton is larger than the standard model result at Q < 600 GeV. However 
the value for the cross section is smaller than that of the standard model. Therefore 
j^SM -g q\q^q |;q j^sm+gr^ loj-gQY values for Q the cross section for the graviton 



'5M 



IS 



is larger than the standard model result but K < K so that also here K 
close to i^'S'A^+Gi? rpj^g i^g^ feature is different from that observed for the LHC in 
Fig. 11a where X'S'M+Gii . This is because the gluon flux does not dominate 

the processes at the Tevatron quite contrary to what is observed for the LHC. For 
the in Fig. 14b and F-distributions in Fig. 14c which are taken at Q = 400 GcV, 
the is close to }{S^^+gr^ Finally we look at the scale variation in Fig. 15a. For 

11 < Ho there is a better improvement when we go from LO to NLO than that observed 
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for the LHC. On the other hand the difference between the graviton mediated process 
and the graviton plus the standard model becomes less. Also this is an indication that 
the gluon flux is less important for the Tevatron than for the LHC. The same feature 
is observed at Q = 400 GeV for the xp distribution in Fig. 15b and the F-distribution 
in Fig. 15c. Because the full X-factor does not differ very much from K^'^ , it may 
appear that it is sufficient to use the SM K-factor in extracting bounds on Ms at the 
Tevatron. Indeed, existing bounds on Ms from the dilepton production process at the 
Tevatron [18] are obtained using a constant 7^-factor of 1.3. These bounds range from 
1.4 to 1 TeV for d=2 to d=6. While the inclusion of the NLO corrections computed 
here may not change the Tevatron bounds significantly, it will certainly stabilise the 
cross-section with respect to scale variations. It is for this reason that the inclusion of 
the NLO corrections becomes important even for the Tevatron. 

Finally we consider the forward backward asymmetry defined by 



This asymmetry is a result of the interferences of the parity violating structure of Z 
exchange with the photon [19] and with the graviton initiated process. The dominant 
contribution to ApB comes from the qq subprocess and hence is important at the 
Tevatron. At LHC for the definition of scattering angle (Eq. 6.2) the ApB will be 
neghgible [19] since the LO contribution vanishes. For the Tevatron the plots are made 
in Fig. 16 for the parameters d = 2 and Ms = 1.5 TeV in the region 20 GeV < Q < 800 
GeV. Here we have plotted Aps ior the SM and the SM+GR at NLO. The effects of 
extra dimensions show up at Q > 300 GeV. In this region the graviton contribution is 
slightly larger than the SM where the latter shows a constant behaviour. 

In summary, we have computed the cross sections da/dQ, d'^a/dQdxp, d^a/dQdY 
and the forward backward asymmetry up to next to leading order. The standard 
model result was already given 25 years ago but now we also included all subprocesses 
due to the graviton in the context of TeV-scale gravity models. Our main conclusion 
is that the NLO QCD corrections are very significant at the LHC because of the large 
incident gluon fiux at the LHC. At the Tevatron where the gluon fiux is small, the 
NLO effects are also quite small. But, significantly, at both colliders inclusion of the 
NLO QCD corrections help stabilise the cross-section with respect to scale variations. 
The extraction of bounds from both colliders will, therefore, require the inclusion of 
these NLO corrections. 
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Appendix A 



Here we present the ajF-distributions which were too long to be pubhshed in section 4. 
The Born contributions (LO) to the coefficient functions are 

27r Hqq[xi, X2, n'p) 



^5M,(0)/ 



N 



<-y»0 I <T^0 



TT Hqg{Xi, X2, H'p) 

8N + xl 



TT Hgg(Xl, X2, H'^p) 

2{N^ - 1) + x^2 



(A.l) 



For the NLO corrections we start with processes with a photon or Z-boson in the 
intermediate state. The corrections are given by [6], [13]- [16] 
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(A.2) 
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l<-*2 



(A.3) 
(A.4) 



We have introduced the following abbreviations 

Hab,12{Xl, X2, 1^%) = Hab{Xi, X2, - Hab{x°i, X2, //p) - Hab{Xi, X^, flj,) 

''rHab{xi, X2, fip) , 

-f^a6,l(^l' A*f) = Hab{xi, Z, /jj,) - Hab{xl, Z, /jj,) , 

Hab,2{^^^2,IJ'F) = Hab{z,X2,IJ,F) - Hab{z,x'^2^fil) , (A.5) 
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The Ca2-i C,})^ and can be obtained from >Ca^, C,\,^ and >Cci by using 1 <-> 2 symmetry. 
The NLO corrections with a graviton in the intermediate state are given by 



Hqq{Xl, X2-I Hp) 



{xi — x\ 



r 

o^ '^ci 



{xl + xl) 



^0 



x\ 2 2 > 

-2 4 + -+ -2 

X\ \ Xi Xi 



ai 



+ j dxi 



(xi - x1){x1 + x^) 



4£ 



+ 



Hgg,l{Xl,X2,l4') 



^2 



+ / dxi / (ia;2 



2 Hqq^i2{Xi,X2,(4p) 
' {Xi — x\) {X2 — X2) {x\ + 0:2) 



Hqq{x^, X2, l^p) 



+6 C(2) - 10 



(3 + 2£(a;?,a;0))log 



+ {C{xlxl)f 



+ 



X 



dxr> 



Hqq(xi,X2, 1^^) 



I "^"^^ / ^"^^ (xi + X2)3(a;i + xl) {X2 + X?) (x? + 
^ (2x^x?a;^ + (x? + x^) (2a;^a;?x^ + Sx2x\xf + Sxf )) 

+— (2xlixl + x°) + 2x1 (xf + + 7xf) 

Xi ^ ^ ' 

+4x2X° (\\xf + 7x?x° - 4xf ) + %x\x% (3xf + - 2xf)) 

,02 



5X 



2^2 



1 f^O 



? + x^) + ( 28x? - 16\ + 44a;^ + 0:2 ( 8 - 16^ 

X2 







X 

X2 



+66a;f + 2x2 (| - 7a;? + 8^ - 7a;0 + j + 2 (Sxf - Ixfx^ 



33 



-7a;?a;f + 8xf) 



+ 1^2 , 



(A.7) 



TT 



4 



dxi 



Hgg{Xl, X2, fJ^'p) 



{xi - xl) 
+ / dxi / dx2 



dxi 



+ x^) 

Hgg^l{Xl, X2, fMp) 



02 

4 ^ + 44 

JU-^ JU-^ JU\ JU 



8 4 \ 
I ) 



{x^-xl){x\ + xl) 

Hgg^l{Xl, X2, /il') 



4£, 



61 



^02 
^2 







8 4 



(xi -a;?)(a;? + a;^) 



-4^ + 4^-- + - 

•^2 -^2 -^S 



(a;? + X^) 



^ + 2 xO) j log (^^) + 



+6 C(2) - 



203 



18 



4, /g^ 

- o log — 

3 V/^F 



35 



Hgg{Xl,X2, l-ip) 



I ^fTfHgg{x°^,xl,nl) 

+0^1^! dx^ J dX2-^^ + X2)3(xi + X0)(X2 + + 0:°) 

^ (^4xlxfx'2 + {xfxf + xfx^^ + 4x2(xf + a:f ) 
\ Uxlixf - xlxl) + xl (Axf + 12xf - Ax\xf) 



+■ 



+0:2 ( - 4.xlxf + IQxfxl + 12a;f ) + Sxfxf + Sx\ 



03^02 



+ — ( 3^2 (-4a;? + Sx^) + X2 ( - 32a;?xf - 16xf + 12xf + 8xf ) 

Xl \ 



34 



( - 32x^x^2 + 8a;f - I6xf) + ^ (l2xfxf + Uxfxf) 



+8x'ixf + 200:^0;^^ + 12x^^"x^2 1 + ^ + I 12 



.02 n2 



4xf 



Xn 



,02 



+xl U'^ - 28xl \ +xi i- 64x°a;^ - 12xf 



Ax 



o3 



24x1 



o2 



Xr) 



+x, ( - 4^ - 4x? - 4:^ - 4x0 ) + ( - 8p - 12p + 40 



+4x^ ( 1 - 1 ) ) - Mx^xf - \Axfxl + 2xf + 2x\ 



X-i X' 



+ 1^2 , 



(A.8) 



(x? + xl) 



.02 



Xa x^ \X2 



Xn 



X 



02\ 



4 D 4" 4 g I 



X2 J 



+ 



+ 



jdx, I 



dxo 



Hqg^l{Xi, X2, h'f) 



(xi - Xi)(a;" + X2) [x2 X2 



^02 

xi 



,^0 
X2 



Hqg{xi, X2, l-lp) 



.0 



16 16 



^ + 8^ ^ ~o r^^i 

Xi \ Xi X\ XI J 



+ 



dx 



jdx, J 
+ / dxi / (ix2 



Hqg,2(Xi, X2, fJ^'^p) 



^x? _ 16 16 

^ (X2 — X2)(X° + X2) xf Xi X? 

H,^g{Xl,X2.JJ^p) 



{xi + a;2)2(a;i + 0:2) (a;2 + x1){x1 + x^) 

— ( — 9t^0™03 _ r, o2 02 I ™2 / _ o 0™0\ 

+X2 ( - 8xf x° - 8x?xf ) + - ( - 2xf xf - 2xfxf) 



+— I - 2xf - \2xXxl - \Ax\xl^ + x\ (164' - 8a;'i' 

X\ \ 



.02^0 



.0^.02 , ,-v.2 



-,.0 _ 0^.0 



35 




36 



Appendix B 

Here we present the ^-distributions which were too long to be pubhshed in section 5. 
The Born contributions (LO) to the coefficient functions are 
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The /Cqj, /Cfoa can be obtained from /C^^, /Cf,^ and /Cc^ by using 1 4^ 2 symmetry. 

The NLO corrections with a graviton in the intermediate state are given by 
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Figure Captions 



Fig. 10a The invariant di-lepton pair mass distribution da/dQ for = 2 TeV and 
d = 3 at = 14 TeV (LHC). Standard model (dotted line), gravity (long 
dashed line), standard model plus gravity (solid line). 

Fig. 10b The double differential cross section d'^a/dQ dxp where xp is the longitudi- 
nal momentum fraction of the di-lepton pair. The plot is given for Mg = 2 TeV 
and d = 3 at = 14 TeV and Q = 700 GeV (LHC). Standard model (dotted 
line), gravity (long dashed line), standard model plus gravity (solid line). 

Fig. 10c The double differential cross section d'^a/dQ dY where Y is the rapidity of 
the di-lepton pair. The plot is given for = 2 TeV and d = 3 at = 14 TeV 
and Q = 700 GcV (LHC). Standard model (dotted line), gravity (long dashed 
line), standard model plus gravity (solid line). 

Fig. 11a The X^-factor for the cross section da/dQ at M, = 2 TeV and d = 3. The 
plot is made for the LHC (v^S* = 14 TeV). Standard model (dotted line), gravity 
(long dashed line), standard model plus gravity (solid line). 

Fig. lib The K^-factor for the double differential cross section d'^a/dQ dxp at Mg = 
2 TeV and d = 3. The plot is made for the LHC with a/S* = 14 TeV and 
Q = 700 GeV. Standard model (dotted line), gravity (long dashed line), standard 
model plus gravity (solid line). 

Fig. 11c The '^-factor for the double differential cross section d^a/dQ dY at Ms = 
2 TcV and d = 3. The plot is made for the LHC with = 14 TeV and 
Q = 700 GeV. Standard model (dotted line), gravity (long dashed line), standard 
model plus gravity (solid line). 

Fig. 12a The ratio Eq. (7.7) for the cross section da/dQ at Mg = 2 TcV and d = 3. 
The plot is made for the LHC (v^ = 14 TeV). Gravity alone in LO (dot-dashed 
line), standard model plus gravity in LO (short dashed line). Gravity alone in 
NLO (long dashed line), standard model plus gravity in NLO (sohd line). 

Fig. 12b The ratio Eq. (7.8) for the cross section d^a/dQ dxp at M, = 2 TeV 
and d = 3. The plot is made for the LHC with ^/S ^ U TeV), Q = 700 GeV 
and Xp = 0. Gravity alone in LO (dot-dashed line), standard model plus gravity 
in LO (short dashed line). Gravity alone in NLO (long dashed line), standard 
model plus gravity in NLO (solid line). 

Fig. 12c The ratio R^ Eq. (7.9) for the cross section d'^a/dQ dY at Ms = 2 TcV 
and d = 3. The plot is made for the LHC with = 14 TeV), Q = 700 GeV 
and Y = 0. Gravity alone in LO (dot-dashed line), standard model plus gravity 
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in LO (short dashed hne), Gravity alone in NLO (long dashed hne), standard 
model plus gravity in NLO (sohd line), 

Fig. 13a The invariant di-lepton pair mass distribution da/dQ for Mg — 2 TeV and 
d = 2 at = 1.96 TeV (TEVATRON). Standard model (dotted line), gravity 
(long dashed line), standard model plus gravity (solid hne). 

Fig. 13b The double differential cross section d'^a/dQ dxp where xp is the longitudi- 
nal momentum fraction of the di-lepton pair. The plot is given for Ms = 2 TeV 
and = 2 at = 1.96 TeV and Q = 400 GeV (TEVATRON). Standard model 
(dotted line) , gravity (long dashed line) , standard model plus gravity (solid line) . 

Fig. 13c The double differential cross section d^a/dQ dY where Y is the rapidity of 
the di-lepton pair. The plot is given for Mg = 2 TeV and d = 2 aX y/S = 1.96 TeV 
and Q = 400 GeV (TEVATRON). Standard model (dotted hne), gravity (long 
dashed line), standard model plus gravity (solid line). 

Fig. 14a The A'^-factor for the cross section da/dQ at Ms = 2 TeV and d = 2. The 
plot is made for the TEVATRON {y/S = 1.96 TeV). Standard model (dotted 
line), gravity (long dashed line), standard model plus gravity (solid line). 

Fig. 14b The A' '^-factor for the double differential cross section d'^a/dQ dxp at Ms = 
2 TeV and d = 2. The plot is made for the TEVATRON with = 1.96 TeV 
and Q — 400 GeV. Standard model (dotted hne), gravity (long dashed hne), 
standard model plus gravity (solid line). 

Fig. 14c The AT^-factor for the double differential cross section d^a/dQ dY at Mg = 
2 TeV and d = 2. The plot is made for the TEVATRON with = 1.96 TeV 
and Q — 400 GeV. Standard model (dotted hne), gravity (long dashed hne), 
standard model plus gravity (solid line). 

Fig. 15a The ratio Eq. (7.7) for the cross section da/dQ at = 2 TeV and 
d=2. The plot is made for the TEVATRON (v^ = 1.96 TeV). Gravity alone 
in LO (dot-dashed line), standard model plus gravity in LO (short dashed line). 
Gravity alone in NLO (long dashed line), standard model plus gravity in NLO 
(solid line). 

Fig. 15b The ratio Eq. (7.8) for the cross section d'^a/dQ dxp at Ms = 2 TeV and 
d = 2. The plot is made for the TEVATRON with y/S = 1.96 TeV, Q = 400 GeV 
and xp = 0- Gravity alone in LO (dot-dashed line), standard model plus gravity 
in LO (short dashed line). Gravity alone in NLO (long dashed hne), standard 
model plus gravity in NLO (sohd line). 
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Fig. 15c The ratio Eq. (7.9) for the cross section (Pa/dQ dY at M, = 2 TeV and 

d = 2. The plot is made for the TEVATRON with 7^ = 1.96 TcV, Q = 400 GeV 
and Y = 0. Gravity alone in LO (dot-dashed line), standard model plus gravity 
in LO (short dashed line), Gravity alone in NLO (long dashed line), standard 
model plus gravity in NLO (solid line). 

Fig. 16 The forward backward asymmetry ApB in NLO for the TEVATRON = 
1.96 TeV. The parameters are Mg = 1.5 TeV and d = 2. Standard model (dotted 
line), standard model plus gravity (solid line). 
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Fig. 10a 
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Fig. 10b 
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Fig. 10c 
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Fig. 11c 
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Fig. 12b 
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Fig. 12c 
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Fig. 13a 
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Fig. 13b 
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Fig. 14b 
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Fig. 14c 
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Fig. 15b 
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Fig. 15c 
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Fig. 16 
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